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Although,  the  problem  of  ablative  and  evaporative  cooling  in  the 
interior  ballistics  had  Been  considered  before,  the  -mechanism  of  the 
cooling  is  not  clearly  understood.  In  order  to  optimize  the  cooling 
effectiveness  a  thorough  theoretical  and  experimental  study  must  be  made. 

This  report  summarizes  the  effort  ;ln  the  period  from  January  1970 
to  April  1971  to  study  the  cooling  mechanism  in  the  interior  ballistics. 

The  basic  aim  is  to  explore  a  better  prediction  of  film  cooling  convective 
coefficient  for  a  simplified  model  which  however  retains  the  basic  charac¬ 
teristics  of  the  interior  ballistics  problem. 

The  experimental  study  was  carried  out  at  Rock  Island  Arsenal  by 
personnel  from  Science  and  Technology  Laboratory  of  the  Rock  Island  Arsenal. 
This  part  includes  design  and  modification  of  cooling  projectile,  conver¬ 
sion  and  interpretation  of  experimental  data.  The  theoretical  study  deals 
with  the  solution  of  gas  dynamics and  heat  transfer  in  a  unsteady* incompres-r 
sible,  two-phase  flow  with  non-isothermal  wall.  Heat  transfer  convective 
coefficients  are  predicted  or  correlated  with  the  experimental  results. 

The  work  reported  herein  was  done  under  contract  number  DAAF-01-C-0380 
with  Rock  Island  Arsenal,  U.S.  Amy  Weapons  Command. 
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NOMENCLATURE 


Dimensional  Nomenclature 


>* 


D 

c 

k 

k„ 


Mj 

P 


Total  cross-sectional  area  of  the  small  holes  cn  the  projectile 
Constant  pressure  specific  heat  of  the  propellant  gas 

r 

Constant  pressure  specific  heat  of  the  liquid 
Constant  volume  specific  heat'  of  the  liquid 

Internal  diameter  of  the  barrel 

Thermal  conductivity  of  the  gas 

Thermal  conductivity  of  the  liquid 
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Gas  temperature 
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Reference  temperature 

Wall  temperature  without  cooling 


V 


Ti&e. 

T*. 

t-t 

P 

v.<rrespondent  time  *.?here  the  p?  vjectile  is.  located 
Reference  time 

Gat'  velocity^  in  radial  direr-ion 
Exit  Telocity  of  the  projectile 
Liquid  velocity'  in  R-direction 
Reference  vleocity  (  U£) 

Velocity-  in  ;  ■'direction 

Gas  velocity  in  Z-direction 

Liquid  velocity  in  z-direction 

Projectile  velocity 

»  y^e  (Rq-R> ,  distance  from  the  wall 

Position  of  gas-liquid  interface 

*  -Y,  pointing  to  the  solid  wall 

corrdiante  along  the  barrel 

Position  of  the  projectile 

Gas  density 

Liquid  density 

Reference  desnity 

Gas  yiscosity 

Liquid  viscosity 

Coordinate  around  the  center  of  the  barrel 

»> 


vl 


if 


ef 

fc 
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ft 
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Non  Dimensional  Nomenclature. 

.  n~^- 
A  n  ,  a  conatant 

a  Total  cross-sectional  area  of  the  snsll  holes  in  the  projectfl 

C*  »  0.5928,  a  constant 

C  «»  0,72 6,  a  cose canc 

dQ  Internal  diameter  or  the  barrel 

E  Eckert  :.u$0er  of  the  gas 

E^  Eckert  mmSe ;  of  the  liquid 

f^  Functions  of  n 

^2  Functions  of  n 

^3  Functions  of  n 


A  function  of  (t,y) 

Heat  transfer  coefficient 

Thickness  of  the  liquid  at  z  »  z 

P 

■  1.82,  a  conatant 
Total  amount  of  the  liquid 
*  -3.65  ,  a  constant 
Nusselt  number 
Pressure  of  the  gas 
Pressure  at  the  hreech 

Pressure  of  the  liquid 
Pjandtl  number 

Heat  generation  in  the  boundary  layer 
Heat  generation  in  the  core  solution 

vil 


s&aaaii 


w 

1  jp 

1 1-  . 

I 

f  r 

F  r 
i 

Coordinate  in  radial  direction 

§  9 

pi* 

k- 

£  c 

1"  >  | 

1 

>*0 

E 

Internal  radius  of  the  barrel 

Reynolds  number 

^  * 

Time 

I  ‘ 

E 

Gas  velocity  in  f'-direction 

S  p 

r  u« 

Liquid  velocity-  in  f'-direction 
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«v 
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M 

U 

■  -  /Re  u 

i 

Sr' 

ik  * 
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Gas  velcotiy  in  z-direction 

g 

1' 

| 

Liquid  velocity  in  z-'direction 

l  - 
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“p 

Velocity  of  the  projectile 

:■■ 
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»  /Re  Cyq-y) 

Thickness  of  the  liquid  film 
coordinate  along  the  barrel 


Position  of  the  projectile 
Gas  density 

Similarity  transformation  invariant  variables 


Gas  temperature 

Gas  temperature  in  the  core 

Liquid  temperature 

Wall  temperature  with  cooling 

Wall  teniperature  without  cooling 

Thermal  diffusivity 
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I,  INTRODUCTION 

It  is  v^ell  known  that  the  present  firing  arms  are.  limited  for  their 
continuous  firing  capability  due  to  the  thermal  expansion  or  cook-off  from 

Jk» 

the  excessive  heating  in  firing.  For  example,  Cohn  Cl)  and  Corner  (2) 

2 

show  that  a  heat  input  of  1,000  Btu/ft  sec.  is  possible.  Many  approaches  have 
been  attempted  to  remove  this  large  amount  of  heat  transfer  such  as  the  change 
of  configuration,  different  propellant,  change  in  firing  frequency,  and 
liquid  cooling.  Without  changing  the  existing  design,  one  of  the  most 
promising  methods  of  cooling  is  to  coat  an  ablative  and  evaporative  material 
on  the  interior  wall  By  the  projectile  as  it  passes  through.  The  thermal 
resistance  of  the  so  generated  vapor  may  seal  off  the  heat  transfer  from  the 
combusted  gas  to  the  wall  .  It  was  shown  by  Cohn  (.13  and  Adams  and  et  al  C.3) 
that  even  smearing  with,  silicone  oil  or  coating  with  teflon  or  beeswax  on  a 
projectile  does  have  some  cooling  effect.  However  the  exact  mechanism 
of  the  insulation  is  not  clearly  understood.  The  present  theoretical  and 
experimental  research  is  motivated  to  study  the  cooling  mechanism  and  to 
optimize  the  effectiveness  of  such  cooling  devices. 

The  experimental  study  used  an  XM140  Aircraft  Automatic  Gun  as  the  model 
to  analyze  the  cooling  effectiveness  of  the  cooling  projectile  which  contains 
water  as  coolant  and  to  obtain  the  boundary  conditions  that  are  needed  in  the 
theoretical  analysis.  The  theoretical  study  analyzes  the  gas  dynamics  of  the 
propellant  gas  to  the  barrel  occurring  behind  the  projectile. 

Convectivs  heat  transfer  coefficient  which  is  important  in  heat 
transfer  calculation  was  deriyed  for  unsteady  compressible  flow  with  and 
without  liquid  cooling  in  a  barrel. 

*  The  number  in  the  bracket  denotes  the  reference  number  listed  in  REFERENCES 
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II  CONCLUSIONS 


Before  details  of  experimentation  and  theoretical  analysis  are  reported, 
important  conclusions  are  first  presented  in  this  section  for  readers  who 
are  interested  in  the  present  results  and  the  application  of  the  solutions 
obtained.  It  is  concluded  that: 

1.  Solid  material  coated  on  the  projectile  is  not  an  effective  coolant 
because  of  the  low-  melting  rate  in  a  time  interval  of  two  to  three 
milliseconds.  Therefore,  liquid  is  adapted  as  the  coolant. 

2.  Water  as  a  coolant  is  proved  to  be  excellent,  since  it  is  inert, 
nontoxic  in  both  liquid  and  vapor  phase,  more  transparent  than  silicone 
oil  in  vapor  phase.  Antifreezer  may  be  added  to  water  to  lower  the 
freezing  point  in  use. 

3.  Proper  geometry  of  the  modified  projectile  that  contains  coolant  is  essential 
for  an  effective  and  uniform  cooling.  Of  the  eight  modified  configurations 

(see  Pig  II-l)  tested,  the  No.  2  type  gave  a  distinctively  effective  and 
uniform  cooling.  These  modifications  are  based  on  the  analysis  given  in 
Appendix  IIA.  Ve  conclude  that  injection  of  coolant  through  side  copper 
band  and  the  front  portion  of  a  projectile  gives  a  better,  even  distribu¬ 
tion  of  coolant  than  other  combinations. 

4.  Experimenting  with  water  as  the  coolant  in  the  No.  2  type  modified  projectile 
the  peak  temperature  on  the  interior  wall  is  reduced  by  as  much  as  40%  in  com¬ 
parison  with  the  standard  round  based  on  the  single  shot  experiment  as 

shown  in  Fig  11-2,  Cooling  is  particularly  effective  near  the  muzzle  end. 

5.  The  analysis,  of  the  gas.  dynamics  of  the  propellant  gas.  behind  the  projectile 
and  the  heat  transfer  through,  the  propellant  and  the  coolant  to  the  wall 
may  be  divided  into  three  regions  -  core  flow,  gas  boundary  layer  flow 

and  liquid  layer  flow. 
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6.  Gas  and  liquid  boundary  layers  are  found  to  be  small  compared  with  the  diameter 
of  the  har*el.  The  thickness  of  the  layer  is  about  one  hundredth,  of  the 
diameter, 

7, "  The  gas  core,  solution  is  solved  provided  that 

the  pressure  response  at  the  breech,  and  the  projectile  position  as  a 
function  of  time  are  known.  The  result  giyes  the  following  dependent 
Variables  in  dimensionless  form  (deviation  given  in  Chapter  V). 

Pressure : 

=tyt>  -  o.(fu6  i  i  -  /( t  (2-D 

Density: 


Temperature : 


(2-2) 


eivt)=  LV«  *  A.*¥#l  -  (ip^-  vu>  )/u  y^}]  (M) 


Gas  velocity: 


W  >  / y-t) 


(2-1*) 


Heat  generation; 


Vip  >  ^  A)  -j>  i‘F  J(f(cy  , ]  (2-5) 

Ty* — LH/H  + 


The  analysis  shows  that  the  pressure  distribution  can  be  predicted  to  3 


6 


15%  accuracy  at  any  location  behind  the  projectile.  Fig  II- 3  shows  the 
comparison  between  theory  and  experiment.  The  aboye  solutions,  are  also 

programmed  in  Appendix  II B  for  general  input  of  projectile  position  versus 
tinse  and  breech  pressure  functions.  The  above  analysis  should 


be  used  0.15  milliseconds  after  the  firing  »  since  the  solution  is 
inaccurate  during  the  strong  propellaut  reaction. 


8.  The  gas  boundary-  layer  flow  without  liquid  cooling  is  first  solved  to 

understand  the  heat  transfer  mechanism.  The  solution  is  obtained  through, 
similarity  transform.  The  analysis  deals  with,  an  unsteady  compressible 
laminar  flow  with  heat  generation  and  non-iso thermal  wall.  For  the 
•  wall  temperature  of  the  type  (see  Chapter  V,  Ea  (5-13)) 


L-  h>  ~  \ t o)]  (2-6) 

the  local  Nusselt  number  which  is  function  of  time  and  position  is 


i.fes- 


ki  _  _  ne'Si  ^  ^ 


(2-7) 


For  the  XML40  Model  or  the  like  f^(o)  =  3.2,f,' (o)  «  -0.12.  It  should 
be  noted  that  Newton’ s  cooling  law  is  defined  here  as 


l-^Tr 


(2-8) 


2. 

-  °r 

where  Tr  * 

and  that  the  time  and  location  in  Eqs  (2*6)  and  (2*7)  naist  be  behind  the 

f  1 

projectile  or  z  <  (t) ,  For  other  type  of  interior  ballistics  -f^CcJ/f^Co) 


in. the  equation  (2-71  ,mav  he  taken  as  a  constant  to  he  determined  from 
experimental  Jnta.  The  limitation  of  the  formula  (2~7)  is  that  the  non- 
isothermal  vail  temperature  must  be  given  by  Eq.  (2-6)  which  may  only 
approximately-  represent  the  actual- response  in  experiment.  Nevertheless 
it  is  recommended  over  the  use  of  steady  incompressible  laminar  formula. 
Examples  of  -using  Eq.  (2-7)  can  be  found  in  Chapter  VI. 

9.  The  heat  transfer  from  the  propellant  to  the  wall  through,  a  liquid  layer 
is  solved  approximately-  in  Chapter  VII.  The  Nusselt  number  which  is  also 
function  of  time  and  distance  along  the  axial  direction  is  given  as 


Nu 


(2-9) 


where  y^  is  the  liquid  film  thickness.  A  method  of  calculating  y^  is  given 
in  Chapter  VII  section  1.  Because  there  is  not  enough  information  about  the 
existence  of  a  liquid  film  behind  the  projectile  and  about  the  percentage 
of  liquid  in  the  modified  projectile  that  may  get  behind  the  projectile  the 
solution  thus  obtained  from  Eq.  (2-9)  can  only  Be  taken  as  preliminary 
one.  Further  study  is  needed  to  verify  many  assumptions  which  are  not 
theoretical,  but  arc  experimental  and  yet  to  be  established.  An 
alternative  heat  transfer  formula  is  also  derived  in  Chapter  VII  section  3 
based  on  the  response  of  the  solid  wall  of  the  barrel  . 

10,  In  general  with  cooling  liquid  there  is  about  30%  of  reduction  in  heat 
flux  from  the  propellant  gas,  to  the  wall  for  the  modified  projectile  of 
No.  2  (see  Fig  11-1)  *  The  mechanism  of  cooling  effect  mainly  comes  from 
the  absorption  of  heat  by  the  liquid  and  therefore  the  reduction  of  temperature 
gradient  near  the  wall. 


Ill  EXPERIMENTAL  STUDIED 
II't-1  Projectile  Modification 

In  the  design  of  the  modified  projectiles  for  experiment,  consid¬ 
erations  were  made  to  provide  an  even  distribution  of  coolant  on  the  interior  wall 
and  a  good  lubrication  effect  for  the  motion  of  projectile.  Then  several 
single  shot  experiments  were  conducted  to  determine  the  best  coolant  exit 
locations  and  hole  diameters  among  the  proposed  modified  models.  The 
hole  diameters  are  calculated  from  a  theoretical  analysis  under  the  cri¬ 
terion  that  the  coolanc  is  completely  squeezed  out  by  the  time  when  the 
projectile  leaves  the  muzzle  end.  In  the  analysis  the  effective  pressure 
to  squeeze  the  coolant  out  of  the  projectile  is  assumed  to  be  one  half  of 
the  bore  pressure  at  that  given  instant  .(See  Appendix  II  A) 

A  capsule  containing  water  is  inserted  into  the  hollow  space  of  the 
projectile  as  shown  in  Fig.  II-l  for  an  initial  test.  One  hole  at  the 
bottom  and  four  to  eight  holes  evenly  spaced  around  the  projectile  are 
drilled.  As  the  projectile  is  fired  the  pressure  on  the  base  of  the 
projectile  will  squeeze  the  water  out  of  the  projectile  through  the  holes 
around  the  projectile.  The  water  then  is  coated  to  the  inner  wall  of  the 
barrel.  The  precise  dimensions  for  machining  are  shown  in  Fig.  II-l. 

Approximately  17  grams  of  water  were  sealed  in  the  projectile  by  a  polyeth¬ 
ylene  capsule,  a  rubber  bag,  and  a  polyethylene  tape  respectively.  The 
experimental  results  show  that  ths  modified  projectiles  NO.  2,  gave  the 
best  result  with  uniform  cooling  and  lower  wall  temperature  (40%  reduction 
in  the  maximum  wall  temperature  from  that  of  the  standard  round.  See 
Comparison  of  Figs.  II-2)  It  was  established  that  the  front  holes  and 
peripheral  hv-^s  located  between  the  copper  band  gave  the  best  result. 

Therefore  No.  2  is  adapted  for  experiments  of  continuous  firing.  The 
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data  obtained  for  single  shot  experiment  are  the  /measurements  of  the 
interior  and  exterior  wall  temperature,  and  the  pressure  response  in  five 
locations  of  the  gun  barrel  .  All  experiments  were  performed  at  the  Rock 
Island  Arsenal  by  personnel  of  the  Science  and  Technology  Laboratory.  The 
data  are  included  in  Appendixes  III  A,B,  and  C. 

1II-2  Instrumentation 

The  Instrumentation  of  the  experiment  was  also  set  up  by 
personnel  of  Science  and  Technology  Laboratory  at  Rock  Island  Arsenal. 

Tb^.  model  used  for  the  experiment  was  XML40  barrel.  Five  pressure  probes, 

fou/c  interior  wall  temperature  probes  and  six  external  thermocouples  were 

used  to  record  data.  Kistler  pizo-electric  high  pressure  transducers  were  used  for 

pressure  measurements.  The  platinum  -platinum  10%  rhodium  Mo-Re  surface 

temperature  probes  were  used  to  measure  internal  wall  temperature.  Both 

probes  have  a  response  time  of  microseconds  which  is  sufficiently  accurate 

for  a  test  interval  of  2  to  3  milliseconds.  The  external  temperatures  were 

measured  by  chromel-Alumel  thermo c  >uples  welded  on  the  outer  barrel  surface.  The 

arrangement  of  probe  locations  is  shown  in  Fig.  III-l. 

III-3  Results  and  Discussion 
HI-31  Recommended  Model 

Eight  proposed  modifications  of  the  projectile  shown  in  Fig.  II-l 
were  nested  on  a  single  shot  basis  in  the  first  two  experiments  and  compared 
with  the  standard  round.  Each  modified  projectile  contains  C.039  lbm 
(17.6  grams)  of  water.  All  eight  modified  projectiles  exhibited 
substantial  cooling  effect  in  the  first  experiment.  Ir.  particular,  the 
peak  internal  surface  temperature  can  be  reduced  as  much  as  40%  from  the 
standard  projectile.  This  can  be  seen  from  figures  ic  Appendix  III.  where 
Internal  surface  temperatures  are  designated  by  S^,  S2,  and  S^.  Among  the 
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eight  experiments  of  modified  projectiles  the  first  three  rounds  tfl,  112 , 
and  j?3  ([specified  in  Pig,  IT-1)  are  judged  to  provide  satisfactory  results. 
For  test  projectiles  #4  through  #8  data  obtained  in  a  short  time  after  th( 
firing  are  satisfactory.  However,  because  of  the  separation  of  the  aluminum 
cap  from  the  projectile  at  approximately  1  millisecond  after  firing,  data 
thus  obtained  are  insufficient  to  indicate  the  effectiveness  of  the  modifi¬ 
cation  Nevertheless,  these  data  provide  a  relative  comparison  among  the 
projectiles  #4  through  if  8.  Based  on  the  results  of  projectiles  ill ,  #2,  and 
//3  projectile  $2  gives  the  best  cooling  effect  through  out  the  barrel. 
Thi3  modification  shows  that  a  better,  uniform  coating  of  liquid  on  the 
surface  of  the  barrel  can  be  achieved  by  squeezing  liquid  out  from  the  front 
portion  and  between  the  copper  rings  of  the.  projectile  as  shown  in  the 
Fig  II- 1  by  P  and  F.  Since  the  cap  of  projectile  separated  from  the  pro¬ 
jectile  for  #5  (which  has  the  same  configuration  as  #2  except  a  larger 
diameter  at  the  bottom  of  the  projectile  denoted  D  in  Fig.  11-1)  this  series 
of  the  test  did  not  allow  us  to  compare  the  effectiveness  of  projectiles 
?2  and  if  5.  Therefore,  the  type  If 2  is  recommended 
III-32  Temperature  Response 

From  Fig.  II-2  we  see  that  with  liquid  cooling  the  peak  internal  wall 
temperatures  are  greatly  reduced.  This  may  attribute  to  the  liquid  sealing 
effect  that  prevents  the  hot.  propnllant  gas  from  contacting  the  wall. 

It  should  be  noted  that  with  a  pressure  higher  than  the  critical  pressure  of 
water  in  the  barrel  the  evaporation  effect  of  liquid  is  almost  negligible. 
Only  when  the  projectile  has 'left  the  muzzle  can  evaporation  be  a  cooling 
'  ffect.  Note  also  that  at  a  given  instant,  say  2.3  millisecond,  Fig.  II-2 
gives  a  high  temperature  for  a  larger  distance  from  the  breech.  This  can 
be  explained  from  the  fact  that  a  larger  projectile  velocity  near  the 


nuzzle  end  contributes  to  a  larger  friction  heating  and  hence  a  larger 
temperature  response. 

Referring  to  Appendix  IIIA,  on  external  wall  temperature  we  observed 
that  ' imperature  response  at  Ks  appears  to  be  the  highest  one.  This  is 
because  the  barrel  thickness  at  K5  position  is  the  thinest  over  the  entire 
length  of  the  barrel.  We  note  also  that  temperature  response  at  is 
high  although  the  thickness  of  ^  is  thick.  However  this  may  be  explained 

by  the  fact  the  core  temperature  inside  the  ^  position  is  hottest  of  the 
core  temperatures.  ,  _ 

III-33  Pressure  lteb{)onse : 

In  the  case  of  liquid  cooling  the  peak  pressure  of  P^,  at  the 
breech 3  was  cut  down  about  5  kpsi,  compared  with,  a  standard  round,  while 
the  response  at  other  positions  remain  the  same.  This  can  be  seen  from 
the  pressure  curves  of  PROOF  ROUND  NO.  3  and  PROJECTILE  NO.  2  tested  on 
June  3,  1970,  Appendix  III  A.  This  is  due  to  a  pressure  release  to  fill 
up  the  hollnw  space  inside  the  modified  projectile  at  the  early  stage. 

However,  this  does  not  seem  to  slow  down  the  exit  velocity.  It  can  be 
assn  from  the  following  fact: 

For  standard  projectile  (without  cooling  PROOF  ROUND  N0.3)  the  pro¬ 
jectile  velocity  at  the  exit,  Wp,  is^;  based  oc  the  core  solution 

i 

W  »  2130.8  ft/sec 
P 

and  for  modified  projectile  (modified  Projectile  No.  2)  its  velocity 

in  ft/millisec.  is, based  on  experimental  data, 

d*  o  '  , 

y  -  — -  1.87696  t  -  0.55565t  +  0.05672t  J 

p  at 

at  the  exit  of  the  barrel  t  «  2.38  millisec. 

W  -  2079  ft /sec. 

Thu  above  evidence  indicates  that  the  exit  velocity  of  the  modified  pro¬ 
jectile  can  reach  the  velocity  attained  by  the  standard  projectile  even  though 
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there  is  a  drop  in  peak  prepare.  This.  indicates,  that  the  liquid  squeezed 
from'  the  lateral  and  front  aidest  of  the  projectile  on  the  wall  of  the 
barrel  does  serve  as  a  lubricant  and  reduces  the  friction  force  between 
the  wall  and  the  projectile.  This  reduction  in  drag  force  provides  a 
higher  acceleration  near  the  mizzle  end  for  the  modified  projectile. 

III-34  Continuous  Firing 

Twenty  rounds  of  the  modified  projectile  number  2  were  tested  in  an 
experiment  of  continuous  firing.  However  only  ten  rounds  were  fired  when 
a  malfunction  of  test  gun  occurred.  The  result  of  the  experiment  was 
given  in  Fig.  III-2  III-3,  III-4,  and  III-5.  Where  Fig.  III-2  gives  the  lo¬ 
cation  of  temperature  probes  and  Fig.  III-3  gives  the  internal  wall  temp¬ 
erature  measurement  at  2.37"  from  the  muzzle  end  without  cooling.*  Figs. 

HI-4  and  III-5  are  external  wall  temperature  measurements  with  and  with¬ 
out  cooling  respectively.  It  should  be  noted  that  there  were  twenty  rounds 
of  standard  projectile  fired  during  the  experiment  which  does  not  allow  a  direct  com¬ 
parison  with  the  data  of  modified  projectiles  where  only  ten  rounds  were  fired. 
Nevertheless ,  we  observed  that  with  and  without  cooling  the  external 
temperature  increases  linearly  with  time  after  the  first  five  rounds. 

Although  based  on  the  single  shot  experiment  the  heat  transfer  behind  the 
projectile  can  be  reduced  by  liquid  cooling,  the  combined  total  cooling 
effect  under  continuous  firing  needs  further  investigation  before  con¬ 
clusions  can  be  reached. 


internal  wall  tempera 


Fig.  ni-fj  External  Vbll  Temperatures  v.s.  Time  for  Continuous 
firing  of  Standard  Projectiles 
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*y  ypWfJLATION  0?  THEORETICAL  ANALXSIS 


IV-'l  General  Consideration 

Consider  that  projectiles  are  continuously-  fired  v*ith  a  certain 
frequency.  The  analysis  of  the  gas  dynamics  Behind  the  projectile  is  then 
basically  a  viscous,  compressible  periodic  flow  with  heat  transfer.  Although 
the  analysis  of  gas  flow  without  liquid  cooling  is  not  contained  in  the 
original  proposal  it  is  included  here  because  with  some  modification  the 
solution  can  be  extended  to  the  flow  with  liquid  cooling. 

The  mathematical  analysis  consists  of  two  essential  parts.  They  are 
the  gas  core  flow  where  the  viscosity  and  conductivity  are  not  important 
and  the  boundary  layer  flow  where  the  viscosity  and  conductivity  plays 
the  major  role.  This  division  is  based  on  the  order  of  magnitude  analysis 
of  an  unsteady  viscous  flow  in  the  following.  From  the  characteristic  of 
the  viscous  diffusion  of  an  unsteady  flow  it  is  known  (for  example  see 
Chapter  V  of  Reference  (4))  that  thickness  of  the  velocity  boundary  layer, 

V*  ___ 

—2- ,  at  a  given  time  under  an  accelerated  flow  is  /  yt  .  Where  £  is 
J.  J  ^ 

the  characteristic  length,  t,  the  characteristic  time,  and  ,  v.  the  kine¬ 
matic  viscosity  cf  the  propellant  gas.  Also  it  is  known  (for  example  see 
Chapter  XII,  p.  270  Effect  of  Prandtl  Number  of  Ref. (43)  that  the  temper¬ 
ature  boundary  layer  is  of  same  order  as  the  velocity  layer  if  the  Prandtl 
nunfeer  of  the  fluid  is  of  the  order  of  magnitude  one.  Consider  that  for 
the  present  problem  £  is  the  radius  of  the  barrel,  0.6  inches^ and  t,  the  duration 

of  the  projectile  in  the  barrel,  2  milliseconds.  Then  the  maximum  boundary 
layer  thickness  is  only  one  hundredth  of  the  radius.  Therefore  99%  of  the  gas 
flow  near  the  center  of  the  barrel  is  not  affected  by  the  viscosity  and  the 
condictivity  of  the  fluid.  Hence  a  core  analysis  may  be  performed  independently 
of  barrel  wall  conditions.  When  the  core  solution  is  obtained  the  boundary 


layer  solution  can  be  solved.  Also  when  the  cooling  liquid  film  is  in- 
eluded, the  analysis  in  the  boundary  layer  becomes  a  rather  complicated 

two-phase  flow.  Fig.  IV-1  illustrates  the  division  of  the  analysis  in 
the  following  sections. 

IV-2  Assumptions 

According  to  the  general  physical  understanding,  the  following  assump¬ 
tions  are  made. 

(1)  The  flow  is  not  too  far  away  from  the  laminar  flow  and  . 

the  Navier-Stokes  Equations  can  be  used'  for  the  present 
analysis . 

(2)  The  gravitational  force  is  negligible. 

(3)  The  flow  is  axially  symmetric  (i.e.,  »  0)  and  there  is  no  circumferential 

velocity  component  (i.e.,  V  «  0). 

(4)  The  liquid  layer  is  incompressible,  and  the  gas  is  compressible. 

(5)  The  flow  is  periodic  in  time. 

(6)  The  wall  is  smooth  but  non-isothermal. 

(7)  Transport  properties  are  constant. 

IV-3  Governing  Equations 

It  will,  be  convenient  to  take  the  cylindrical  coordinate  system  with 
axis  being  the  center  line  of  t;he  barrel.  Let  (U,V,W,)  be  the  velocity 
components  in  the  direction  of  (K,+,Z)  respectively.  ThenU,W,T,P  are 
functions  of  (t  ,  R,  Z)  (see  Fig.  IV-1).  For  the  gas  stream(without 
subscript) : 

(A)  Continuity  equation 

(H) 
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(B)  Momentum  equations 


Rr-directicn, 


nW^sk  -i) = -%~a* 


(Ji-2) 


where  p  is  the  viscosity  of  the  gas,  and 


l-v  Z  +•« 


where  2  and  £  are  unit  vectors 


V  *f  -  -i-  h _ (  C> t r \  ,  ^VT 


Z-direction, 


e  (14-3) 


(C)  Energy  equation 


-  ?§• 


o.  Tir^  r  . 

+>u  SL  +  6lU 


(MO 


where 


4  -  *WEsW«|*Jl .  ts  »  Ci|_(RU)  jgf  (k-a 


k  is  the  thermal  conductivity  of  the  gas,  is  the  heat  generation 


due  to  combustion. 


■»^*  J&  w  '»*^5v' jfc>* 
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(D)  Equation  of  state 


£  ~-f  K  TCc 


(h-6) 


where  C  is  a  correction  factor  depending  on  the  range  of  pressure 


under  consideration.  For  a  pressure  below  37  kpsi  C£  is  0.726  approxi¬ 


mately. 

For  the  liquid  layer  (with,  subscript  l ) : 
(A)  Continuity  equatlcn 


r.  cr(RcjjO-»-  I^W* 


(la-7) 


(B)  Momentum  equations 
R-direction, 


- -$f +/$$&%)  * ^  .  *% 


fc? 


(a-8) 


where  and  p^  are  density  and  viscosity  of  the  liquid  layer  respectively 
Z-direction, 


f<  1 « ^  )  -  -  i  +*  i 


’<>*'  1 


(U-9) 


(C)  Energy  equation 


v»  -  *<  L-kfeom)  -  %) 


( h— 10) 


T„  and  kff  are  temperature  and  thermal  conductivity  respectively. 


Here  we  have  nine  unknowns  (i.e. ,  liquid  U  ,  W  ,  P  ,T  ,  and  gas  U,W,P, 

X.  X*  X.  X 


T,p)  with  nine  equations, 


•  fw-  -  ir  - 
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IV-4  Boundary  and  matching  conditions 

(1)  Boundary  Conditona 

At  wall  (R  »  R  ) 
o 

“t/tl-t  ,  RCy  2.)  -O 
T) t  ?  Tl^  > 


(Ml) 

(1.-12) 


At  breech 


2  -  .) 


W -  W  '^  R .  afu.)) 

HI  -=  o 

T  --  TU,R..V'») 

At  the  projectile  at  any  instance  Z  «  Z^(t) 

P*  Pp  ^ 

“vvy  -'vwf 

■u  --  0 


( 1.-13) 


(h-iU) 


T  -  Tit  ,  R  f  2plO> 

(2)  Matching  conditions  (the  subscript  i  denoting  the  gas-liquid 
interface) 


(A)  Kinematic  condition 

The  continuity  of  velocity  components 

(B)  Dynamic  condition 

The  continuity  of  normal  stress 

T».  -  V!t  lit). «  P.  -  »/*  )A 

The  contiituity  of  tangential  stress 


(h-l?) 


(a- 16) 


;  2.-17) 
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(.C)  Thermodynamic  conditions 
The  conservation  of  energy 


( It— 18) 

The  continuity  of  temperature 

Ti.  •  Ta  <^W> 

where  Ah  is  the  heat  of  evaporation.  However  Ah  is  negligible 
for  the  operating  pressure  near  or  above  the  critical  pressure 
of  the  fluid.  For  water  the  critical  pressure  is  3.21  kpsi. 
IV-5  Nondimensionalizatlon 


We  make  the  above  variables  dimensionless  by  letting 

.HL  a 


l_rr 


Wj 

ur 

4,  -  it-. 

*  Wr 

e  -  3i_ 

01  Tr 

1' 

u'r 

\>  * 

(Vo*' 

-if 

*- 

2 

-t  ^  3„ 

f '  i 

(li-20) 
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where  U^,  py,  L,  tf ,  Ty,  are  dimensional  characteristic  quantities 

defined  as  follows  with  typical  values  in  the  parentheses. 

L  --characteristic  length  of  the  barrel  (3.478  ft) 

-projectile  velocity  at  exit.  (2130.8  ft/sec) 

t  -characteristic  time  —  (1.6322  millisec.) 

U 

r 

p  -  density  of  the  propellant  gas  right  behind  the  projectile  at 

r  the  exit  (5.504  lb  /ft3) 

ra  ,,  2 


U 


characteristic  temperature 


-  (2173.9  °R) 


r  ‘  R 

D  -  diameter  of  the  barrel  (1.2  inch) 

P  -  characteristic  pressure  p^U^.2  (5.388  kpsi) 

Mq  -  molecular  weight  of  the  propellant  gas  (23.805) 


R  -  Gas  constant  (- 
Ah  -  latent  heat 


1545.33  ,,  ftlbf-. 

m - M-9142 

o  m  .R 


with  known  transport  properties 
lb 

_c-  _ m 

u  *  2.8x10  ftsec  (see  TABLE  IVA) 

1  *  »•*«  JK 

«  0.071x10  ftsec 

K  *  0.04  Btu/hr.ft  °R  (1500  °F)  (gaa) 

-  0.35  Btu/hr. ft°R  (500°F) 

cp  -  0.435  Btu/lb°R  (gas) 

*  0.35  Btu/lb°R  (gas) 

Define  the  dimensionless  parameter: 

U  p  L 

Reynolds  nunfcer  R  ■  — r-  — — 
e  V 

CM  C  P 

Prandtl  nutrber  P  »  — ,  P  * - ^rr - 

r  K  ’  r(  K 

II  2  ,  U  2 

Eckert  nunfcer  E  ■  —  ,  E„  *  — ~= - 

Vr  -  11  Cp»Tr 

We  have  then  the  governing  equations,  in ''dimensionless  form  as 
(1)  For  the  gas  stream 

(A)  Continuity  equation 

o  (i-2i) 

(B)  Momentum  equations 


Wit**?  <r/» 

w  ^r,,,.,^,^.,  ■  „  .  . .  .... * 


1 
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TABLE  IVA  Viscosity  of  equilibrium  air 
(  from  ARS  Journal,  Aug.  1961,  p.1152) 


•""■  m  (ntillipoisc) - - - , 

T(*K)  0.1  Mm  1  atm  lOat-.n  100  Mm  5000  atm 


3000 

0.87 

O.SC 

1000 

1.13 

1.12 

5000 

1.37 

1.3G 

0000 

1.G7 

1.C1 

7000 

1.97 

1.91 

8000 

2.20 

2.17 

0.86  0.85  0.85 
1.09  1.03  1.07 
1.31  1.32  1.29 
1.58  1.51  1.51 
1.82  1.78  1.71 
2.11  2.03  1.97 
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CC)  Energy  equation 


e  lH  *  tt57  *wlt  ]  “  p&  ^T  +  E  L5t +a^  *“Hi  u'a) 
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(D)  Equation  of  state 
P  "  ^®cc 

(2)  For  the  liquid  layer 

(A)  Continuity  equation 

T  7?™  ^  * 

(B)  Momentum  equations 
Redirection 


Z- direction 

(dW*  ^ 


(C)  Energy  Equation 

/  ,  .,  i>&> 


(3)  Boundary  and  matching  conditions 

At  wall  (r  *  fQ) 

Uj  ■»  UG*  s  c 

Qfc  * 

At  breech  (z  =  0) 

P  *  pb(t,r,o) 


w  *  0 


u  -  0 


0  -  0(t,r,f') 


(U-25) 


(It- 26) 


(U-27) 


._L_ikT  ^.1  ,,  oa, 

Vi-c  ^  i>r  •>*>■)  G%  " R«  /V  L*  <x*  1  U-28; 


(U-30) 


(V31) 
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At  projectile  (z  *  z^) 
p  -  pp(t,r,  y 
v  -  wp(t,r,zp) 
u  -  0 
0-0 


At  interface 


ii 


V  Vii 


(  k  -*•  -  (  b  4-  —  £M.\ 

vrt  Xk;  ST-1;  '  (P  r« 
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Obviously  the  above  problem  is  difficult  to  solve 
with  the  division  of  flow  region  given  in  IV-1  we 
to  solve  the  problem  in  the  following  sections. 


(4-32) 

(4-33) 

(U-3b) 

( *i-35) 

However 
may  proceed 


V  ANALYSIS  OF  CORE  SOLUTION 


V~1  Core  Problem 


From  the  general  consideration  in  Section  IV-1  we  found  that  the  boundary 
layer  thickness  is  very  small  compared  with  the  diameter  of  the  barrel.  The 
propellant  gas  flow  near  the  center  of  the  barrel  can  be  solved  by  neglecting 
the  boundary  layer.  Accordingly,  the  analysis  becomes  that  of  one-dimensional 
unsteady  compressible  flow  with  the  governing  equations  as  follows: 


(pm)  =.  o 

n  »vv,: 

>U5  .  £b 
+^7rr)  -  ~r~ 


v 


>ctv^e 


(5-D 

(5-2) 

(5-3) 


(5-U) 


where  p,w,p  and  0  are  unknown  core  density,  velocity, pressure, and 
temperature ,  q^  is  the  heat  generation  of  the  propellant  gas. 

The  above  governing  equations  are  solved  with  the  following 
boundary  and  initial  conditions 

at  z  *  z^(.Q)  w  *  0  p  =  {^(t)  (breech,  pressure) 


at  z  *  z  (t) 
P 


V?  »  Wp* 


dz 

_ E. 

dt 


(5-5) 


where  z^(t)  is  the  projectile  position  at  the  given  instant.  For  the 
present  case  z^Co)  ■  0.04642. 


At  Zp  *  0.268  (where  data  is  available)  p  =  2.21 
where  p  is  determined  in  Appendix  VA- 


(5»6) 


at  t  *  0  w  =  0  (5*7) 

The  projectile  position  as  function  of  time  is  obtained  from  the 
experiment  through  the  method  of  least  squares  as  (standard  round) 

zp  *  0.04642  +  0.92578  t2  -  0.11475  t3  -  0.00919t4  (5*8) 

Thus  we  have  the  projectile  velocity 

dz«  o  -K 

1.85156  -  0.34425t  -  0.03676tJ  (5-9) 

p  dt 

Equations  (5_8)  (5-9)  are  plotted  in  Fig.  V-l  and  Fig  V-2. 

V-2  Method  of  Solution 

We  assume  that  the  density  is  a  fvnction  of  time  only.  This  assumption 

leads  to  a  sufficient  condition  of  Lagrange’s  assumption  which  has  been  shown 

by  Heinev  (5)  to  be  an  adequate  one  for  the  interior  ballistic  problem.  Under 

this  assumption  the  pressure  and  temperature  are  still  functions  of  the 

position,  z,  and  time,  t.  From  Eq.  (5“1)  we  have 

3w  1  3p  ,  . 

3z  *  "  p  3t  ’  p  ~ 

Integrating  with  respect  to  z 

Noting  that  at  z  *  0,  w  =  0,  \o  that  fn(t)  =•  0,  and  that  at  z  *  zp,  w  «  wp,  we 


(5-8) 


(5-9) 


w 

-** 

z 

p 


(5-10) 


dt  o  _  1  dp 
z  p  dt 

P 

Integrating  Eq  (5 -10a)  with  respect  to  t 


(5-10a) 


f  1 

P(t)  -  — 
P 


where  C’  is  integration  constant  and  is  determined  as  C*  =  0.5928  by  using 
the  condition  in  Eq  (5-6) « 


Fig.  Projectile  Displacement  and  Acceleration  v.  s.  Time 
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Then 


p(t) 


0.5928 


zp(t) 


(5-11) 


To  find  the  pressure  distribution  we  substitute  Eqs.  (5~10a)and(5-ll)  into 


(5-2)  and  integrate  to  obtain  ,  , 

C-  •;  (/-*>» 


P(z,t)  »  t>b(t)  - 


z  2  2 

P 


C5-12) 


after  satisfying  the  boundary  condition  (5-5  )  .  The  temperature  profile 
follows  readily  by  substituting  (5-10a)  ,(5-ll)»  and  (5-12)  into  Eq.  (5-4). 
That  is 

-Ife 


...  2  2 
^2 _  zp  ~2f 


.5928 


2  z 


(5-13) 
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From  energy  equation  (5-3)  we  obtained  the  heat  generation  of  the  propellant 
as 


q  «  - 
tcr«. 


C*(z-Z2  (o)) 
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zp(°)  C' 


(5-14) 


where  the  dot  denotes  the  differentiation  with  respect  to  time  (i.e. , 


dC  ) 


dt 


(  )).  Therefore  we  have  solved  the  core  solution  of  the  interior 


ballistic  problem  in  terms  of  the  projectile  motion,  Zp(t).,  breech  pressure. 


Pb(t),  for  IMR  propellant.  The  solution  is  computerized  in  Appendix  JIB.  The 


propellant  gas  velocity,  w,  behind  the  projectile,  the  temperature,  T  *  QT^, 


the 


density,  p,  pressure  ,  p,  and  heat  generation  q  are  all  plotted  in  figures 


V-2  V-3,  V-4,  XI-3,  and  V-5. 

V-3  Result  and  Discussion 

Although  many  works  on  core  solution  of  the  interior  ballistics  are 
available  such  as  Spurk  (6 j,  Heiney  (5)  Love  and  Pidduck  (7),  Vottis  (8), 
and  Carriere  (9j  the  present  solution  has  an  advantage  that  the  burning  rate  of 
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the  propellant  and  the  friction  between  the  projectile  and  the  barrel  are 
not  assumed,  instead,  they  are  all  absorbed  in  'he  specified  boundary 
conditions  of  projectile  motion  and  the  breech  pressure.  Both  boundary 
conditions  come  directly  from  experimental  measurements.  Hence  the  error 
due  to  assumptions  of  burning  rate  and  friction  effects  are  eliminated.  Also 
the  present  solution  has  an  advantage  of  simplicity  in  analytical  form  over 
the  solution  obtained  from  theory  of  characteristics. 

In  the  theory  of  characteristics  the  solution  must  be  obtained  from 
series  of  numerical  evaluation  of  characteristics.  Furthermore  when 
heat  generation  is  included  in  the  gas  flow  the  characteristic  ajlution  can 
become  very  complicated.  Therefore  cue  present  solution  althov.gL  ..ess  ver¬ 
satile,  is  more  desirable  for  prediction  of  convective  heat  ti-.i.sfer  through, 
the  gas  boundary  layer. 

-In  Fig.  V-2  we  plot  both  projectile  velocity  and  fluid  velocity 
behind  the  projectile.  The  dashed  line  is  the  experimental  data  obtained 
in  Fig.  V-l  while  the  solid  lines  show  the  theoretical  gas  velocity.  -It 
is  interesting/  to  note  that  the  gas  velocity  near  the.  breech  end  first 
increases  and  then  decreases  but  at  a  larger  distance  the  gas  velocity  decreases 
monotonically  with  tim.-*- 

Fig,  V-3  shows  the  predicted  core  temperature  and  measured  wall  temperature 
with  and  without  cooling.  The  difference  in  the  core  temperature  and  wall 
temperature  at  a  given  time  drives  the  heat  flux  from  the  core  to  the 
barrel.  We  note  that  at  any  given  time  the  core  temperature  behind  the 
projectile  always  decreases  from  the  breech  to  the  muzzle  end.  However  the 
wall  temperature  at  a  given  time  with  and  without  cooling  increases  frot-  the 
breech  to  the  muzzle  end.  More  will  be  discussed  in  Chapter  VI  when  boundary 


39 


layer  analysis  is  given.  Fig.  V-4  gives  the  density  variation  with  respect 
to  time.  The  dashed  at  the  small  time  was  computed  from  the  core 
solution.  However  in  this  period  because  propellant  combustion  is  strong 
it  is  felt  that  the  core  solution  should  not  be  taken  too  seriously  for 
its- validity.  Fig.  V-5  is  the  plot  of  heat  generation  in  the  core.  It  is 
obvious  that  generation  is  the  largest  near  the  breech  end.  About  80  per¬ 
cent  of  total  heat  generation  is  confined  in  the  ten  percent  of  the  barrel 
length  near  the  breech. 

With  the  core  solution  obtained  the  heat  transfer  and  ga3  dynamics  of 
boundary  layer  flow  behind  the  projectile  may  be  analyzed. 
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VI  BOUNDARY  LAYER  SOLUTION  WITHOUT  COOLING 
VI-1  Governing  Equation  and  Similarity  Transform 

As  already  discussed  in  Chapter  IV  there  are  two  flow  regions  in 

the  barrel.  One  is  the  core  stream  in  the  central  portion  of  the  barrel 
and  the  other  is  the  Boundary-  layer  near  the  wall.  This  boundary-  layer 
flow  will  be  analyzed  here.  According  to  the  boundary  layer  theory 
(see  Ref.  C4))  the  pressure  in  the  boundary  layer  is  still  that  of  core 
solution.  The  other  dependent  variables  such  as  density,  temperature, 
velocity,  and  neat  generation  must  now1  be  a  function  of  an  .additional 
variable  r  which  does  not  appear  in  the  core  analysis.  The  core  solution 
how  becomes  the  outer  boundary  condition  of  the  boundary  layer  while 
the  wall  condition  is  the  other  condition  of  the  flow.  Only  when  this  flow 
is  solved  can  the  heat  transfer  between  the  propellant  gas  and  the  barrel 
be  predicted.  This  problem  is  much  more  complicated  to  analyze  in  comparison 
with  steady  boundary  layer  flow  because  there  are  three  independent  variables 
r,  t,z  and  both  outer  and  wall  condition  are  non  isothermal  and  unsteady. 

In  dealing  with  this  complicated  problem,  we  note  that  Reynolds  number 
prUr'"  S 

in  the  flow  is  Re  ■ - »  1.45  x  10  Cor,  in  the  conventional  definitions 

*  PrUrRQ  P  ? 

Re  =  — - - s  2.9  x  10  )  which  is  very  large.  Consequently,  the  boundary 

layer  thickness,  in  order  of  — ~  ,  is  very  small.  Therefore,  for  the 

vRe 

boundary  layer  flow  we  shall  define  the  new  variables 

y  *  /Re  (r  -  f),  distance  from  the  internal rwall  of  the  barrel,  and 

,u  «*>^e  u,  velocity  in  y-direction  (6-1) 

to  stretch  the  small  quantities  (r  -  r)  and  u  so  that  each,  term  in  the 

governing  equations  in  these  new  variables  jLa  of  order  of  unity.  Upon 

neglecting  the  terms  of  small  order  we  reduced  the  Eqs.  (4-21)  through 

(4-24)  in  terms  of  new  variables,  y  and  u  ,  t  .  _  , 

*  *  -  ,  of  gas  boundary  layer  as 
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Continuity  Equation 
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Momentum  Equations 

In  y  -  direction 
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In  z-  -  direction 
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(6-6) 


Equation  of" State 

p.  »  P0  0.726  (see  Eq.  (5-12)  for  p) 

Considering  the  solution  form  obtained  in  the  core  solution  we  see 
that  the  independent  variable  z  can  be  separated  out  in  the  boundary 
layer  flow  by  setting 

w(t,y,z)  *  H(fc,y)z  ( (6-7) 

Furthermore,  the  heat  generation  qn  .  in  Eq.  (6  5).  is  assumed  to  possess 

o«L« 


the  similar  form  of  q  in, Eq.  (5-3),  that  is 
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And  it  is  noted  that  the  Eckert  number 
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in  the  present  case.  The  term  E(-|* —  )  is  thus  negligible  compared  with 

3  1 

other  terms  in  Eq.  (6“5)  which  are  of  order  of  unity. 
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Substituting  Eqs.  (6-6),  (6-7),  and  (6-8)  in  to  Eqs.  (6-2)  to  (6-5)  we 
have  the  resulting  governing  equations 


It-  +  ly  “ :p  *  p1Ln  ° 


$T  +“If  +  h2>  cc’  ?2+|t 
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These  unsteady  compressible  boundary  layer  equations  are  solved  with  initial  and 
houndary  conditions  to  be  described  in  VI-2.  The  method  of  similarity  trans¬ 
formation  will  be  used  to  solve  the  above,  problems.  The  method  of  similarity 
transformation  can  be_  divided  into  two  kinds,  i.e..,  via  separation  of  variables  • 
and  une-parameter  group  theory.  The  excellent  references  for  the 

latter  method  can  be  found  in  the  books  by  Hansen  (IQ)  and  Ames  (11)  and 
the  paper  by  Morgan  (12) . 

z 

From  the  present  experimental  result  we  found  that  in  equation 

at 

z  P 

(6-10)  could  be  approximated  fay  a  form  of  «  mtn  as  shown  in  Fig-  VI-1, 
where  both  m  andtji  are  constants.  This  mahel  it  possible  to  apply  one- 
parameter  group  theory  of  similarity  transform  to  reduce  our  partial 
differential  equations  into  a  set  of  ordinary  differential  equations  which, 
can  be  solved  numerically.  In  order  to  satisfy  the  invariant  requirement 
of  the  method  we  have  to  choose  the  similarity  variables  as  follows 
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where 


_  .  . ^  n  are  new  variables  in  another  group 
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A-  *  -  ,  n  is  the  index  of  mtn 

fj_*  f2»  and  are  arbitrary  functions  of  n. 


Substituting  Eq.  (6-12)  into  Eqs  (6-9),  (6-10),  and  (6-11)  we  have  a  set 


of  oridinary  differential  equations 


(f2  + 


n)  £[  +  fx  f ;  +  (f3  +  n+2)fx  -  0 


(6-13) 


Cflf2  +  *fl>  f3  +  (f3-1)flf3  “  f3  “  «** 


(6-1 h) 


u'l  fl2+  pr  *2*1  '  £Il£i  -  0  (6-15) 

8 

where  the  primes  denote  the  differentiation  with  respect  to  n.  The  detailed 

derivation  is  presented  in  Appendix  VIA.  From  Fig  VI-1  we  see  that 

,2 

a  x  2  n  *  _3  65 

- 5^-  /z  «  mtn  can  be  approximated  by  1.82t  ,  i.e.,  m  *  1.82,  n  *-3.65. 

dt^  p 

These  values  give  A  *  1.325.  The  above  differential  equations  then  can  be 
rewritten  for  XM140  model  as 


(f2  -  1.325  n)f^  +  fjf^  +  (f3  -  1.65)  fx  -  0 


£6-13a) 


(r,f,  -  1.325  nfx)fj  +  (f3  -  1)  M,-f"  *  1.08 
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VI- 2  Initial  and  Boundary  Conditions 

It  is  seen  that  five  conditions  that  combined  initial  and  boundary  conditions 
are  needed  to  solve  the  problem. 

At  the  wall,  y*0  ,  or  n  *  0,  the  density  related  function  can  be  approximated 
from  the  experimental  data(see  Appendix  VI  B)  as 
fj/i)  -  3,2 

and  the  no-slip  and  impermeable  condition  give 


f,<o)  -  0 


(6-16) 


for  u  «  0 


f3(o)  »  0 


for  w  *  0 
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From  tha  core  stream,  y  -*•  «•,  q.-*  «•',  nq  gradients  in  y  direction  can  exist, 


I;  i. 

|  & 


$ 

r  *** 


f  I  L 


If  ii 


f  ’  i  ' 

A; 


£{  C-l  -  0 


0°)  »  Q 


for  f-  -  0 

,  3w-  _ 

for  r — «  0 


C6-171 


Vi-3  Method  of  Solutions 

The  IBM  360  CSMP( Continuing  System  Modeling  Program)  computer  Jinguage 

was  employed  to  provide  the  numerical  solutions.  As  CSMP  program  handles 

initial  value  problem,  only  condition  (6-16)  at  0  can  j,e  msed. 

Therefore,  we  must  guess  the  values  of  f^Co)  and  £3  Co)  such  that  f^(»)  *  0 

and  f^C®)  *  0  at  n  +  «>  *  After  several  attempts  were  made  we  found  that 

f|(o)  *  -0.12  and  f^Co)  *»  2.024  provided  us  satisfactory  results.  The 

computerized  print-plots  are  presented  in  table  and  figures  in  Appendix  V3C. 

\le  examine  the  convergence  of  solutions  in  the  following  manner.  From 

Eq  (6-7)  we  have  p  *■  t^  f^Cn)*  t  *4^f^Cn).  Now  as  n  -*•  00  the  density 

C* 

in  the  boundary  layer, p,  should  approach  Pcore  »  -  ,  where  C'  *  0.5928  and 

1  65  P 

Zp  «  0.589t  *  approximated  from  Fig.  IV— 1.  Thus 

.  .  x-1.65  .  ,  *  0.5928  .-1.65  -1.65  • 

p(t,»)«t  fx(»)  -  o^gg—  t  =t 

that  is,  f^(®)=  1  at  n  +  ®.  This  condition  is  satisfied  by  the  numerical 
solution  as  shown  in  the  first  figure  in  Appendix  VIC. 

From  Eq.  (6-13a) ,  as  n  +  ®,  f^  +  0,  then  f^  *  1.65.  From  the 

computerized  print-plots  it  is  seen  that  f^  converges  to  -2.426  where 
f j  0  satisfying  =  0,  and  f^  converges  to  1.65,  as  n  goes  larger. 

Finally,  Eqs.  (6-12)  and  (6-7)  can  be  rewritten  as 
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With,  these  expressions  and  the  table  in  Appendix  VIC  we  are  able  to  evaluate 
the  values  of  p ,‘  u,  and  w-  for  any  given  value  of  n. 

Now  we  examine  the  ^validity  of  our  prior  Boundary  layer  approximation,  in 
which,  we  asstime  the  boundary  layer  is  thin  compared  with  barrel  diameter. 
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To  evaluate  the  maximum  boundary  layer  thickness  that  occurs  at  the  maximum 
time  t  »  1.409  .dimensionless,  (the  time  the  projectile  exits)  we-  observe 
from  .plot  in  Appendix  VIC  that  the  velocity  approaches  the  core  solution 
fj  -*■  1.65  at  n  -*■  4.  That  is  in  y  dimensionless  coordinate 

y  «  t1,325  n  »  C1.409)1,32x  4  «  6.3 
The  boundary  layer  thickness  is  estimated  to  be 
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It  is  readily  seen  that  the  ratio  of  boundary  layer  thickness  to  the  barrel 
radius  is  *|g-  which  is  about  one  percent  of  the  barrel  radius.  Therefore, 
we  conclude  that  our  prior  boundary  layer  approximation  is  valid. 

VI-4  Heat  Transfer 
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Now  we  should  consider  heat  flux  from  the  gas  to  the  cold  wall  and  the 
heat  transfer  coefficient  in  the  temperature  field. 
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From  the  equation  of  state  and  the  expressions  for  p  and  p  in  Eqs(5-12) 


and  C6“19) »  respectively,  we  obtain 

fl  «  P  »  t-1*65  *1  -  r_  ,  Aft--3*65 

0  0«726p  0.726.  fph  ~  1,08t 
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Here  we  should  remark,  that  the  houndary-  condition  f^Co)  *  3.2  in  Eq(6-16) 
was  approximated  one.  For  the  detailed  evaluation  is  referred  to  Appendix  VIB. 
In  this  approximation  we  did  not  take  into  account  of  the  short  time  interval 
in  which  the  temperature  changed  from  room  temperature  to  the  approximated  mean 
value,  750°F.  Therefore,  the  temperature  profile  obtained  from  Eq.  (6-22)  and 
heat  transfer  formulas  must  exclude  that  short-  time  interval. 

Under  the  above  restriction  the  local  heat  transfer  is  given  by  Fourier's 


;  1 
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R  -  R 


(6-23) 


Recall  that 


y  ■  J^Re  (~~  -  ^  )  *  (vo-y) 

Defining  Hu  *  ,  Nusselt  number,  from  the  modified  Newton's  cooling 

1®*  q  «  hTr  we  obtain  the  dimensionless  form  of.  Eq(6-23)  as. 


Nu  ■  l'Re 


(6-24) 


The  expression  for  r—  can  be  derived  from  Eq.  (6-22).  Upon  substituting 

y  *y-0 

we  obtain 


c.iis- 
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Hhctre  fj(p)  -  7*1.12,  fjCol  »  3.2  are  obtained  preyioualy  fjrpm  similarity 
solution.  Thus,  we  haye.  for  laminar  unsteady  compressible  flow-  in  a  pipe 
with,  a  standard  projectile  tHe  local  beat  transfer  coefficient  as 

(6-26) 

Since  the  last  term  in  the  equation  (6—2S)  usually  ie  small  compared 
with.  pb  term  for  approximate  evaluation  the  Nusselt  number  may  be  taken  as 


N*  —  0.0 f 4 /A  -fc  fyfe) 

where 

P.  -m  dimensionless  breech  pressure 

*  %  - 


«r  "  the  exit  velocity 


(6-27) 


Pr  -  the  reference  density  right  behind  the  projectile  at  the  exit. 


VI-5  Discussion 

From  the  f2  print-plot  in  Appendix  VIC  we  note  that  f2(F2  in  computer 
program  notation)  keeps  increasing  till  a  vlaue  of  n  »  1.4  is  reached.  Then 
it  starts  decreasing  and  finally  converges  to  a  value  of  f2  »  -2.426.  As  the 
Vertical  component  of  gas  velocity  u  relates  to  f2  by  Eq(6-20)  this  means,  that 
there  are  two  flows  coming  together  from  the  wail  and  the  core,  respectively, 
and  then  joining  at  a  region  where  n.  »  1.4.  Ke  note  that  the  true  value  of 
the  velocity  u  is  very  small  since  p  relates  to  u  by  is  ■-u/f?£~  .  in  which 
Be  generally  is  very  large.  Therefore,  the  radial  velocity  u  does  not 

affect  significantly  the  one  dimensional  core  flew  outside  the  boundary 
layer. 
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In  the.  following  an  example  and  a  comparison  o|  the  present  result  with, 
incompressible  steady  state  flat  plate  boundary-  layer  solution  may  illustrate 
the  procedure  for  calculation.  Consider  the' heat  transfer  calculation 
‘Kith,  present  formula  at  the  time  t  *  1.2  &.e.  t  »  1.9585  milliseconds) 
and  the  position  z  *  0.5075  Ci.e.  Z  ■  1.7655  ft  from  breech)  we  have  for  the 
Standard  projectile  in  XM  140  model 


?r“5.S0  4 
L  =  3.478  -ft 

-0.962.  Ufc  tM.») 

Tr  *  2 1 74-  CK. 


v,-- i.JO.g 

-  2.8  X  l  0~r 

4t  *«t. 

tv  -  K  04 


Re  ^  K  -  oc,4 


Thus  from  Eq.  (6-26)  we  have 


N„  =  %sS.a2«/03 


the  heat  flux  is  then 

2 

q  ■  hTr  -  403  Btu/ft  sec. 

Although  there  is  no  similarity  for  the  physical  phenomenon  between  the 
present  solution  and  steady  laminar  flow.  If  however  the  laminar  solution 
is  applied  instantaneously  at  a  local  position  some  comparison  may  be  made. 

Ve  thus  now  assume  that  at  the  time  t  *  1.2  and  position  z  *  0.5076  the  steady 
incompressible  flow  solution  (see  p.  285  of  Ref.  4)  applies  the  above  example. 


That  is 


Nn  (  Nj.r) 


(6-28) 


where 


*U  «  1375  ft/sec  (from  core  solution) 
Z  ■  0,5076  x  L 
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N  p  £jsndtl  number  =  0.75 
pX  * 


Nu  «-—«  0.7237  X  104 
z  1c 

heat  transfer  is 

«  ■  fcCIcore  -  ’‘wall1  " 


°r  x  10  -  k  'C3250  -  520)  -  1340 - |£S— 

°-5076  1  fc2aec. 


The  above  comparison  shows  the  steady  state  heat  transfer  over  estimates 
heat  .transfer  at  t  *  1.2  and  z  *  0.507.  Therefore  the  steady  flow  formula 
is  not  applicable  to  the  present  problem. 
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YU  THEORETICAL  ANALYSIS  OR  TWO  PHASE  FLOW 
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In  the  anlysis  of  liquid  cooling  effect  many  difficulties  and  uncer¬ 
tainties  appear.  For  example,  ft  is  not  known  in  what  manner  the 
cooling  liquid  is  squeezed  out  from  the  modified  projectile;  how  much 
of  the  liquid  will  get  behind  the  projectile;  and  -t-feat  whether  or  not  tlie 
liquid  will  stay  in  a  film  layer.  Also  the  extent  of  the  influence 
of  grooves  on  the  cooling  effect  is  not  known.  However  an  attempt  is  made 
to  obtain  some  solutions  even  they  are  solved  under  crude  assumptions.  Therefore, 
the  analysis  herein  should  be  viewed  as  preliminary*  result  and  further 
improvement  certainly  is  necessary. 


VII-1  Determination  of  liquid  Film  Thickness 

Ir.  the  following  analysis  we  shall  assume  for  lack  of  experimental 
evidence  that  the  liquid  squeezed  out  from  the  projectile  will  form  a  film 
behind  the  projectile.  Then  the  total  amount  of  the  liquid  behind  the 
projectile,  (t) ,  from  the  breech  to  the  projectile  at  any  given  time 


t  is 


MjlV)  -  r(R.;-R^d2 

» n 


(7-3) 


where  z  *  axial  length 


»  liquid  density  *  const. 


inside  diameter  of  the  barrel 


Rq-R^  *  thickness  of  liquid  film  on  the  wall 


radius  .of  interface 


Eq,  (7-1)  can  he  written  in  dimensionless  form  as 
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(7-2) 
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Where 


.  .ii. 

{*  n 


v  -  V*.  n 

»o  -  —r~  'a 


Ra' 


_  \>A 


dfc  =• 


^  V4 


*  =  A 


Although,  the  exact  amount  of  liquid,  ^Ct),  at  a  given  time  is.  not  known 
a  crude  estintate- of  n^Ct)  is  given  in  Appendix  VIIA  from  which  ye  have 


™‘it)  -  *  r^  ■**' 


■t  ^  K*e*\ 


(7- 2a) 


where  pCt')  is  the  pressure  at  the  base  of  the  projectile. 

The  liquid  layer  thickness  Cro“^)  may  be  approximated  by  a 


polynomial'  form. 


k-r*-  fy-t) 


(7-3) 


Where  ffl(t)  »  function  of  t,  and 

W  and"a2  are  instants  to  be  determined  as  follows. 


At  the  breech,  z  *  0,  Ci^-r,)  »  0,  so  that  ao  *  Q.  Eq.  (7-3)  becomes 


r4-rA-  - Ca,v*A»V] 


(7-W 


Substituting  Eq  C7-A)  into  Eq(7-2)  and  integrating  we  obtain 

■*, m--  (tni.  iK  L-rV-T^rl 


t*7-S> 


S3 


At  the  position  immediately  hehind  the  projectile,  z  *  z?,  *  h^Ct). 

VJfcfere  b^Ctl  is  the  liquid  thickness  at  the  Base  o t  the  projectile  and  is 
derived'  in  Appendix  VITA.  ?tom  Eq,  (7~4I  we  have 

V>K  4*  t  \  Vj>  *  3  ( 7.6) 


Then  and  are  determined  from  Eqs  (7-5)  and  (7-6)'  as 


CVx  >. 


-  v  Lf; 

^  Vk/i 

6  VvMfeJ 

^xr<ic4^j: 
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(7-7) 


"(k  Vj>v 


^t'n«(o'V«k^rp 


Substituting  Eq.  47-7)  into  Eq  (7-4)  and  noting  that  y±  »  /  Re  Cv»0-r±)  we  have 

P  •— 


Once  mt  (Eq.7-2a),  h^  (Eq.7-6) ,  and  Zp  (Eq.5-8)  are  calculated  the  liquid 
film  thickness,  y-» ,  at  a  given  time  and  a  given  position  along  z-axis  can 
be  evaluated  from  Eq. (7-8) . 

A  alternative  simple  estimation  of  liquid  film  thickness  is  also  presented 
as  follows; 

water  filled  in  the  projectile,  Mq  *  0.039  lB^. 

3 

water  density,  ■  62.4  lbm/ft 
diameter  of  the  barrel,  D0  *  1.2  inches 
length  of  the  barrel,  L  **  3.473  ft 

Consider  that  all  water  is  completely  squeezed  out  and  uniformly  coated  on 
the  wall  of  the  barrel  during  the  time  interval,  t  «  0  ~  2.3  m.s.  Then 
the  liquid  film  thickness,  RQ-R^,  ia 


tiol^ 
(7-8) 


**-**  =  '  TTvn  J^s  =  ^,iix‘°  »*■ 

l* 


The  ratio  of  R  rR  to  barrel  radius  K  is 

•  •  •  °  is  .  o 


=•0.0  U^S* 


whicbL  is  one  percent  of  the  barrel  radius. 


V1I-2  Two— Fhase  Gas  and  Liquid  Layer  Flow 

There  are  three  zones  in  the  flow,  free  stream  (core  solution)  around 
the  barrel  center,  liquid  layer  near  the  wall,  and  gas  boundary  layer  in 
between.  Among  these  regions  the  core  solution  it'  already-  obtained  previously. 

VII-21  Gas  Boundary-  Layers  with  Liquid  Cooling 

The  liquid  layer  from  the  above  estimation  is  very  thin  and  the  flow 
velocity  in  it  is  small  compared  with  that  in  gas  layer.  Therefore,  the 
exi stance  of  the  liquid  layer  does  not  affect  the  gas  velocity  boundary  layer 
presented  in  Chapter  VI.  For  the  gas  temperature  layer,  since  the  cooler  liquid 
film  presents  near  the  wall,  it  needs  some  modification.  We  thus  assume  that 
the  mass  flux  pu  in  the  y  direction  remains  that  of  the  gas  boundary  layer 
without  cooling. 

Substituting  0  *  ^  into  Eq.  (fr^ll)  we  have  the  energy  equation 

=  tO  ,7  . 

f “‘i  Rs?f-  (7'?) 

Making  the  substitution  for  p  and  u  from  Eq.  (6~12)  and  setting 
f^  *  fjfg.  Here  cart  be  approximated  by  the  following 


(7-10) 


»  3?  1  j^P 

We ‘have  t~  fA  *  — .  — j-  ,  A  -  1.325 


**  * 


living  by  •jj^—  ,  i  t  becomes 


i'* 

li 

j  i 

Li 


t'V  s-LilM.L|  $* 

+  1®.  "  pv  ei^^) 

<J 


Integrating  once 


b-f-^  ~  |  4v’^ 

^A 


is  *V  !  ^4j  +  fu.vi 

— —  -=  £*  3  a 

*3 


where  f(tjz)  i3  a  function  of  t  and  z. 


(7-11) 


Integrating  Eq  Cl • 11)  once  more 


A  *V t  1.  +f  ft>*>  , 

-  I  e  *j  +-f-^-M 


As  y  ■*  “>,  0  *  0 


f, where  0  is  given  in  Eq.  (5-13) 


Then  the  temperature  distirbuticn  across  the  gas  thermal  boundary  layer  is 
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(7-13) 


The  Interface  temperature  will  be  determined  by  Eq  (7-37)  later.  Here 
y^(t.z)  is  the  position  of  the  interface,  given  in  Eq.  (7-8). 

V1I-22  Liquid  Velocity  Field 

Since  the  velocities  and  in  liquid  layer  are  small,  we  consider 
viscosity  and  pressure  terms  only.  The  governing  equation  is,  therefore, 
from  Fq.  (4-27) in  y  coordinate  variable 


££  m  Jh_ 

1  M  r}* 
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Integrating  with  respect  to  y 

-  *  t  »  -  J.IV> 
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At  the  gas-liquid  interface  the  shear  force  should  be  matched,  i.e. 


?kvM  M, 

\  A‘  *  3 


Then  g-^Cz)  is  determined  as 


(7-16) 
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Integrating  Eq^(7-i5)  once  .^ore 


Wt  -  - 

>^i'LD  5*  ^ 


"•  *>i  1  ■  £b  } 

M] +  ^v> 


It  is  readily  seen  that  2)  *  0  for  y  ■  0,  v.  -  Q« 
As  the  liquid  layer  is  very  thin  we  tsay  approximate 


r3!i‘  c'-> 


which  is  obtained  Chapter  VI. 

Since  Eq.  (6-7)  gives  w  «  2  H(t,y),  we  have 


1  ui 


iH 


“A-  I 


-2.  Jai* 


4 co)  ^  *.o*4*yt 


where  we  have  substituted  H  ^  t  and  f^(o)  *  2.024. 

Upon  substituting  Eq.  (7-19)  back  to  Eq.(7-17)  and  noting  that 
-J.tr 


-s.tr 

Tr 


(7-17) 


(7-18) 


(7-19) 


(7-20) 


from  the  core  solution  we  obtain  a  final  form  for  w^  in  liquid  layer 
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yXI-23  Liquid  Temperature  Field 

Kith,  the  velocity  profile  known  we  nay  approximately  solve  .the  liquid 
energy  equation  by  assuming  that  the  dissipation  is  negligible. 

We  have  the  enrgy-  equation  as 


»±- A 

at  *3  a  “r 

where  ^  -  ------ 

nwn 


(7-22) 


(-7-23) 


i=  T 


Since  the  equation  is  still  complicated  we  solve  it  by  an  approximation 
method  in  which  the  temperature  profile  is  assumed  in- view  of  2q.  (5-13)  as 


0*  *  C0(  +  C,te.‘ 

where  c  and  are  to  be  determined  with  known  conditions. 

O  4 

We  knew  that  at  z  >  0,  0  ■  B(z  »  0)  because  there  will  be  ho  liquid 

*■  <ctf4 

cooling  or  gas  boundary  layer.  Thus  we  have  C  (t,y)  -  8  Ci"0). 

°  tCt-» 

Eq.  (7-24)  becomes 
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If  c^(t,y)  is  further  approximated  by 
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tempe.rature  prpfile  Kith,  the  gas  profile.  Tp  dp  this,  t*e  us.e  the  matching 
condition  of  heat  flux  at  interface  y  t*  y^,  That  is 
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(7-36) 


Substitute  Eqs.  (7-11)  and  (7-13)  for  the  left  hand  side  of  Eq.  (7-36)  and 
Eq..  (7-32)  evaluated  at  interface  for  the  right  hand  side  of  Eq.  (7-30)  we  have 


ife.  k 


ii 


'  V" 


B  •+ 

PUO* 


ec  —(.lax* 

<i  \  i  ^  ft.  :>  V1  \(  p  ***  ^ 

•  -g-  g-t  *- *  « JAV  21 _ ±^l_ 


(7-37) 


«  ►.  .  -*•***  I  ■» 


F»t 


.  i  r  *•  -  ^ 

lt  is.  suggested  that  this  equation  should  be  evaluated  in  the  further  study. 
VH-3  Discussion 

First  it  should  be  remarked  that  the  solution  obtained  is  an  approximated 
one.  Further  improvement  and  study  is  certainly  needed.  However  some 
conclusion  based  on  the  present  result  may  be  drawn.  Let  us  consider  a 
calculation  of  heat  transfer  with  liquid  cooling  by  Eq.  (7-35)  for 
XM140  Model  at  t  «  1.2  and  z  *  0.507. 
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(7-35) 


Here  y^  is  the  dimensionless  liquid  thickness.  To  estimate  this  thickness 

we  assume  that  three  quarters  of  liquid  in  the  projectile  may  reach  behind 
the  projectile  and  is  uniformly  coated  on  the  barrel  surface.  Thus  from 
Eq.  (7-8) 
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(Cgj*'  0^)  in  Eq,  C7“35)  is  the  dimensionless  temperature  difference  between 
the  gas.  liquid.  -  interface  and  the  wall.  To  estimate  this  we.  assume  that  the 
gas-liquid  interface  temperature  is  approximately  equal  to  that  of  the 
wall  temperature  of  the  standard  round  without  cooling.  From  Fig  IX-2  we  find 
that  at  z  *  0.507  Ci.e.,  where  S£  probe  is  located}  at  t  ■  1.2  (i.e. ,  1.9586 
milliseconds)  the  temperature  difference  is 
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For  — -  ,  the  time  rate  change  of  the  wall  temperature,  we  have  from 

Fig  11-2  at  the  same  time  and  position  as 
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Q  from  Eq.  (7-23)  gives 
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Note  that  Reynolds  number  and  the  reference  quantities  are  the  same  as 

given  in  Chapter  IV.  Now  we  may  calculate  the  local  Nusselt  number  from  Eq.  (7—35) 

as 


Nv*.-  -  4.63  x/o2 


To  calculate  heat  transfer  we  have  from  the  definition  of  modified  Newton's  Law 
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V*  281  Wtc  % 


Recall  that  we  had  q  ■»  403  Rtx  *  sec  for  without  cooling  .  This  means 
that  there  is  approximately-  30%  in  reduction  of  heat  transfer  to  the  wall 
at  t  «  1.2  and  z  »  0,507  when  modified  projectile  was-  -used.  Although 
the  ahoye  calculation  gives  approximately-  the  right,  order  of  reduction  in 
heat  transfer  it  is  felt  that  there  are  many  assumptions  involved  that  may 
not  be  easily  verified. 

An  alternative  short  way  to  estimate  the  heat  flux  and  its  coefficient 
is  presented  as  follows.  From  the  experiment  we  learned  that  the  heat  flow 
does  not  reach,  the  external  wall  of  the  Barrel  at  t  =»  2,3  m.s.  after  the 
projectile  was  fired.  We  may,  therefore,  assume  that  within  this  time  interval 


where 


the  heat  transfer  may  be  represented  by  the  heat  flow-  to  a  semi-infinite  solid 

-*/ 


body  exposed  to  a  time-dependent  temperature  on  the  surface  as  t  >  0,  See  Fig 
50H  The  internal  barrel  surface  is  simulated  as  the  surface  of  the  semi-infinite 
body.  Then  at  a  given  .  "■  .tion  of  z  we  have  the  following  governing  equation 
on  solid  wall  side,  in  dimensional 
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temperature  profile  in  the  solid  wall 


Tq  «  room  temperature  at  t  »  0 


time  variable 
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Ys  *  distance  from  the  surface  pointing  to  the  solid  wall 
a  *  the  thermal  diffusitivity  of  the  solid. 

Let 
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where 


t 

P 


the  correspondent  time  where  the  projectile  is  located,  in  the 
barrel  after  the  firing 


We  have  Eq.  (7-38) 
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with  the  initial  condition 
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s  w  *  *  wow 

Y  ®  T  a  0 
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The  solution  fo;r  T  fs,  readily  obtained  from  the  book  by  Carslav  and  Jaeger 

Cl3)p.62, 
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where  A  is  a  dummy  variable.  The  internal  surface  temperature  Tw(A,Z)  is  a 
function  of  time  and  position  along  Z-axis  and  can  be  approximated  from  experi¬ 
mental  data.  Once  Tw  is  known  the  integration  in  Eq(7.44)  can  be  performed. 
However,  the  irregularity  of  the  present  experimental  curves.  Fig  VII-1,  did  not 
allow  us  to  provide  an  appropriate  equation  for  5.  at  the  time  when  this  report 
was  written.  A  further  effort  is  certainly  needed  to  investigate  the  Eq.(7.44). 
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Then  the  local  heat  flux  to  the  solid  wall  is 


i’-K  45-1 


(7-liS) 
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where 


-Y8  «y  *  distance  from  the  surface  pointing  to  the  gas  flow 
Fron  Newton's  cooling  law  we  have 


t  s  **  ^ ~  T.)  -  ^  llTur*~To )  -Ttj] 
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Kith.  Eq.  (7-45)  we  have 


(7-1*7) 
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Then  the  heat  coefficient  is 


Vc  — 
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(7-48) 
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Substituting  for 


3T 

3Y 


from  Eq.  (3-36)  we  obtain 


i  i 


L 


'tl  -! _ I*  T«U.t) 

s-"r*  1  U'- 


Thus  the  local  Kusselt  number  becomes 
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(7.49) 
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where 


(7-50) 


ka  *  Conductivity  of  solid  wall 
k  *  Conductivity  of  gas  boundary  layer 

$w(tjz)  =  internal  surface  temperature  of  the  barrel  which  is 
given  in  Fig  VII- 1 

Due  tc  the  irregularity  of  the  experimental  curves  it  is  very  difficult  to 
provide  an  approximated  curve  to  represent  Tv(t* ,z)  properly.  However 
numerically  integration  of  Eq(7-50)  may  be  performed  to  find  the  local 
Nusslet  number,  it  should  be  noted  that  for  this  Nusselt  number  the  heat 
transfer  calculation  should  follow  the  conventional  Newton's  cooling  law  given 
In  E9.  (7-46). 
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(a)  Without  cooling 


(b)  with  cooling 

Note:  (1)  T  «  Tv~To,  where  T^  *  absolute  wall  temperature,  TQ  * 
room_teroperature. 

(?)  t*  =  t-tn,  where  t  *  any  time  from  0  to  2.3  millisec., 
tp  *  time  interval  for  the  projectile  passes  a  given 
position  after  firiig,^  t  £  tp. 

(3)  bashed  lines  indicate  t>2.3  millisec.,  which  does  not 
apply  to  the  present  problem. 

Fig.  VII-1  Internal  Wal#l  temperatures,  fw  ,  With  and  Without  0  onling 
v.s.  lime,  t  . 
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yill.  RECOJWENDATI.ON 

The  study  of  the  la3t  fifteen  months  under  the  contract  recommends 
the  following  continuing  studies: 

(A)  The  experimental  results  obtained  up  to  now  show  that  the  Modified 
projectile  with  coolant  dees  have  substantial  cooling  effect.  However, 
much  more  data  is  needed  in  order  to  understand  the  mechanism  of  cooling 
and  to  calculate  heat  transfer.  An  improvement  in  measuring  the  interior 
surface  temperature  is  needed.  In  particular,  at  the  breech  end  a  high 
temperature  probe  with  fast  response  should  be  installed.  This  measurement 
is  important  since  it  is  one  of  the  boundary  conditions  needed  in  the  heat 
transfer  analysis.  .  A  measurement  of  the  core  temperature  behind  the 
projectile  can  .be  made  by  implanting  a  thermocouple  in  the  projectile 
flush  with  its  base.  The  thermocouple  wire  protected  in  plastic  tube 

is  leu  chrough  the  gun  barrel  to  the  muzzle  end  and  is  then  connected 
to  the  recording  device.  To  the  present  investigator's  knowledge  the 
measurement  of  core  temperature  following  a  projectile  ha3  not  been  done 
before . 

(B)  Movie  pictures  may  be  used  to  record  the  spreading  of  the  coolant  when 
it  comes  out  of  the  muzzle  end. 

(C)  Hot  wire  measurements  of  velocity  and  Schlieren  pictures  can  be  used  to 
study  the  gas  dynamic  behavior  of  the  mixing  of  the  coolant  and  the  pro¬ 
pellant  gas  at  the  exit. 

(D)  A  series  of  continuous  firing  should  be  continued  at  the  Rock  Island 
Arsenal  to  determine  the  cooling  effect  of  the  modified  projectile  No.  2, 

That  is  to  find  for  the  modified  projectile  the  maximum  interior  and 
exterior  wall  temperature,  maximum  heat  transfer,  maximum  rounds  cf  continuous 
firing  within  the  limits  of  yield  stress  and  cook-off. 


(E)  From  the  measured  data,  the  friction  between  the  modified  projectile 
and  the  gun  barrel,  a  parameter  heretofore  unavailable  can' be  calculated. 
Hence,  the  efficiency  of  lubrication  by  the  coolant  may  be  determined. 

(F)  The  mathematical  analysis  obtained  under  the  present  contract  is 
preliminary  in  nature.  Further  study  is  certainly  needed  in  order  to 
achieve  a  better  prediction.  Refinement  of  the  core  solution  can  be  achieved 
by  considering  that  the  density  is  both  time  and  spatially  dependent.  The 
present  laminar  solution  of  unsteady,  compressible  flow  should  be  extended 

to  the  turbulent  region  so  that  the  effect  of  the  mixing  of  the 
coolant  and  the  propellant  gas  can  be  included. 

(G)  A  suitable  integral  method  can  be  developed  for  the  turbulent,  unsteady 
compressible  gun  barrel  flow.  This  is  not  the  conventional  K&rm&n  - 
Pohlhausen  method  since  in  the  present  case  there  are  three  independent  var¬ 
iables  namely,  the  axial  variable, z,  the  radial  variable,  r,  and  the  time,  t. 
The  derivation  can  be  made  specifically  for  the  interior  ballistics  problem. 

(H)  Study  of  the  flow  in  front  of  the  projectile  has  little  effect  on 

heat  transfer  between  the  propellant  gas  and  the  gun  barrel.  However,  in  case 
of  continuous  firing  gas  dynamics  in  front  of  a  projectile  may  .appreciably 
affect  the  amount  of  heat  transfer.  These  gas  dynamic  effects  include  the 
shock  formation  and  the  propagation  of  expansion  waves. 

(I)  Eq.  (7-37)  should  be  computed  to  compare  with  the  assumed  experimental 
one. 

(J)  Some  form  of  study  on  the  mixing  of  the  cooling  liquid  and  the  pro¬ 
pellant  gas  should  be  initiated  since  it  is  important  for  the  calculation  of 
heat  transfer. 

(K)  The  effect  of  rifling  groove  must  be  considered  in  order  that  a  better 
prediction  of  heat  transfer  coefficient  can  be  achieved. 
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APPENDIX  IIA 

DETERMINATION  OF  MODIFIED  PROJECTILE 


'nsider  a  total  amount  of  water  **  Q.G39  Hun  In  a -modified  projectile 
to  be  squeezed  out  in  a  time  interval  of  2.3  milliseconds. 

Based  on  Fig. '  Ii*'3  we  choose  a  mean  pressure  difference  »  Ap ,  between  both 
ends  of  the  small  hole  to  be  Ap  =  10.  KPSi  approximately. 

Substituting  into  the  formula 


7=/-f 

We  have 


2g 


,  where  V  is  the  water  velocity  at  the  exit  of  the 
-■small-  -hole- 


-  /  10x10^x144 
62,4 


x  2x32.2  -  1220  ft/sec. 


Then  the  .total  area  needed  for  the  small  holes  is 


2.*  3  X 


-4,  . 

JD  fp- 


If  the  nuniber  of  the  small  holes  is  8,  then  their  diameter  is 


D,  0<?|?  “ 
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APPENDIX  IIIB  EXPERIMENTAL  DATA  (June  19,  1970) 


NOTEi  See  Pig.  III-l  for  Arrangement  of  Probes 
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APPENDIX  yA. 

-  DETERMINATION  of  CONSTANTS  FOR  Eqs.  (5-6)  AND  (5-11) 

Pressure,  see  Fig  II-O  and  Fig  IUrl 

P2  -  10.2  KPSi,  at  z  n  0.26804 

?e  ■  3  KPSI ,  at  z  -1,  i.e.-,  at  exit- 

Temperature,  estimated  from  Fig.  H-2  and  Ref.  C3)  Fig.  9 

Tsl  -  2100°F,  at  z  -  0.26804 

Tge  -.1200°?,  at  z  -  1 

Then  dimensionless  density 
_  '«%“*!*  “ 
p  .  s~  .  -  i  a 

pr  feTsl  -2-2l 
From  Eq.  (5-11)  C’  is  determined  as 
C*  *  pz  «  2.21  x  0.25 SO 4  -  0.5928 
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APPENDIX  VIA 

DERIVATION  OF  EqS.  (6-13)  THROUGH  (6rl5)  EY  $l«ILARITy  TRANSFORMATION 


Let  new  variables  (One  Parameter  Group  Theory) 

ll  --  a-  $ 


'V  'ti 

m  -  a.  h 


*%» 

U  -  a  u 


(6A-1) 


where  a  is  a  parameter,  a^,  a^f  y^  Y3  are  constants .  In  terns  of 
the  above  new  variables  Eqs  (b-9)  through  (G-ll)  are  written  as  follows 


--Vi* 


f  ^  a  -S*  ^ 


pH 


-*s}’  .  ,  ^  a.i^-'IxOrv 

--  ft  Wc'i  H.  ft 

<1 


J»-u-£,  ^>r  »  |sF t-*  ***Jl.£l 

t,«4-  -  7. Wf )  -  *  c  T 


(6A-I4) 


In  order  to  have  a  conformal  invariance  the  index  of  parameter  a  in  each 
term  of  the  equation  oust  be  equal  .  That  is 


*1-"*,  *-  <i>  •*  -V"** 


(6k-S) 


r)  %  ~  ~  “  k\— ’ll  —*'(>*  X  -'tv—  *~iz  -  -v\J,  a  S.J 


(6A-6) 


^ V.  "’ti >  — ~i x.  -  ^  <)v 


(6A-?) 
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(6A-2) 

(&-3> 
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The. constant*  in  the  ahoye  equations,  can  be.  determined  as  follows 


(6A-8) 


frost' Eq.  (6A-33  we  obtain 


“of;  A  -1 


(6A-9) 


from  Eq*  C6A-7) 


13“  *« 


(6A-10) 


from  Eq.  C6A-6) 


—  =  >.A-w 
from  Eqs  (6A-9)  and  (6A-11) 

A  -  *  -  -|  A  - 


■n-l 


( 6A-11) 


(6A-12) 


Now  we  can  define  two  new  invariant  variables  5  and  n  as  follows: 


t. 


-JL. 


.  ,^<=2/ai  ^A 


and'  let  the  dependent  yariable&  be 

rv  ,  a.HA, 

f  *  \  *  4%tt)  ~  t  •ft1'*) 
U  -  f  ^  "f  ty) 

H  *  T*’  4i  l1?) 


(6A-13) 

(6A-lli) 


(6A-15) 
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where. fir  f2,  and  ate- functions  of  rL# 
By  the  chain  rule  we  have 

it  -  c'z  '-»* 


(6A-16) 


of  transformed  ordinary  differential  equation 

(4*-^)+/  +f,-k' +  ('+»•*  =-0 

*+rV’+  fv4v4/-+rv  ~  o 


(6-13) 

I  ' 

(6-HO 

(6-15) 


where  the  parameter  a’s  cancel  one  another. 
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APPENDIX  VIB 


DETERMINATION  Qp  ^(n)  AT  vQ  FOE  EQ.  (6-16) 

From  the  transformation  Eq  (6-12) 

„  ,  -  1.65 

f  ■  P.t 

At  the  well,  y  *  0  ,  n  ■  0,  by  equation  of  state 
c  ,  .  P  ,1.65 

flCo)  "  “vfe726t 


( 6B-1) 


(oB-2) 


The  pressure,  obtained  froa  core  solution  Eq.  (5-12)  can  be  readily  fed  in. 
The  wall  temperature  in  Eq.  (09^2)18  measured  froa  Fig  IX.-2  at  positions 
z  •  0.26804,  0.50763,  and  0.69451,  for  various  tiae,  t,.  Then  f^(o)  at 
different  tiae  are  determined  by  a  coaputer  program.  Hie  resulted  output 
for  fjXo)  (we  used  c  in  the.  place  of  f^(o)  in  the  coaputer  program)  is 
tabulated  in  the  following  coaputer  program.  Since,  f^(o)  should  be  a 
constant  in- equations  (6-13)  through  (6-15),  hence  the  approximated  value 
of  f1(o)  which  is  taken  as  average  of  various  f^o)  is  computed.  This 
resulted  in  f^(o)  -3.2. 
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Appendix  yu.A 

DERIVATION  OF  s^Ctl  FOR  Eq.  (7-2)  AND  h^Ct).  FOR  EQ.  (7-6). 

The  liquid  film  thickness  right  behind  the  projectile  is  h£  and  total  liquid 
flow  rate  squeezed  from  the  small  holes  on  the  projectile  i$ 


(7A-1) 


where 


«  density 

A-  »  total  cross-section  area  of  the  small  holes 

* 

Ap  *  pressure  difference  between  two  ends  of  the  small  holes 


In  dimensionless,  Eq(7A-l)is  written  as 


(7A-2) 


i--p.  24 


:  t  |  — 

id  & 


where  (See  also  Nomenclature) 

*  _  Ma  -  A  * 

***  ^LzTTt 


$l5l. 

3- 


Therefore  the  total  amount  of  the  liquid  from  the  breech,  to  the  base  of 

the  projectile  which  is  at  z  -  zp(t)  is  ^  c\  J  Lflx'i  ^  ^  ^  <l.ec*\ 
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If  Ap  m  |>  fa.  assumed  the  difference  fn  the  pleasure  between  the  base  of 
the  projectile  and  exit  hole  of  the  projectile.  >fe  h$ye 


•'w^K)  = 


f  t  p— ..— 

*  Jm 

A  *  ft, 


<L\ 


•fc  £.  vw'.thui) 


To  estimate  the  rate  of  coating  on  the  wall  of  the  barrell  we  lote  that 


(7A-3) 


where 


d  - 


jr-  -  diameter  of  the  barrel 


ftp  -  velocity  of  the  projectile 


For  conservation,  Eq  (7A-2)  and  dq  (7A-3)  must  be  equal.  Then 


C\ 


^cU:f  i  fv 


(VA-h) 


If  we  assume  Ap  »  |  ,  p  Is  obtained  frts  the  core  solution  Eq  (5~12) ,  Eq.  (D-4) 


is  rewritten  as 


h  •  -- 


lit.: 


TTdl 


J  a 


(7A-5) 
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